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Abstract. Using the ergodicity principle for the expectation values of several types 
of observables, we investigate the thermalization process in isolated fermionic systems. 
These are described by the two-body random ensemble, which is a paradigmatic model 
to study quantum chaos and specially the dynamical transition from integrability to 
chaos. By means of exact diagonalizations we analyze the relevance of the eigenstate 
thermalization hypothesis as well as the influence of other factors, like the energy and 
structure of the initial state, or the dimension of the Hilbert space. We also obtain 
analytical expressions linking the degree of thermalization for a given observable with 
the so-called number of principal components for transition strengths originated at 
a given energy, with the dimensions of the whole Hilbert space and microcanonical 
energy shell, and with the correlations generated by the observable. As the strength 
of the residual interaction is increased an order-to-chaos transition takes place, and we 
show that the onset of Wigner spectral fluctuations, which is the standard signature 
of chaos, is not sufBcient to guarantee thermalization in finite systems. When all 
the signatures of chaos are fulfilled, including the quasi complete delocalization of 
eigenfunctions, the eigenstate thermalization hypothesis is the mechanism responsible 
for the thermalization of certain types of observables, such as (linear combinations of) 
occupancies and strength function operators. Our results also suggest that fully chaotic 
systems will thermalize relative to most observables in the thermodynamic limit. 

Keywords: Connections between chaos and statistical physics, Matrix models, 
Quantum chaos 
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1. Introduction 

The success of thermodynamics is based on the fact that the state of a macroscopic 
system behaves as it were independent of the microscopic details of the system, and 
thus it can be determined from a few universal laws. In classical isolated systems 
this thermodynamic universality can be explained satisfactorily if one assumes that 
the dynamics is ergodic and mixing. These properties define classical deterministic 
chaos [I] and lead to the equiprobability of equal- volume regions of the available phase 
space, which is the starting point of the microcanonical ensemble formulation. This 
implies that, on average, the state of the system is independent of the initial conditions 
as well as independent of the measurement time. It is thus generally accepted that 
the system equilibrates into a thermal state, described by the microcanonical ensemble, 
when its dynamical regime is ergodic. On the other hand thermalization is expected 
to be inhibited for integrable or quasi integrable systems. In these systems (almost) 
every trajectory in phase space lies in one of the so-called invariant tori, which foliate 
the phase space. Therefore, being restricted the trajectories to these structures, the 
equiprobability equal- volume regions of the available phase space does not hold and the 
system will not thermalize, at least in the sense explained above. 

In recent years the study of the equilibration and thermalization mechanisms 
in isolated quantum systems has attracted a great interest partly because the 
non-equilibrium dynamics, after an external perturbation is applied, has become 
experimentally accessible for ultra-cold quantum gases [2], and electrons in solids . 
The technology makes it possible to induce sharp changes in the parameters controlling 
the system and then observe the subsequent time evolution, which is essentially unitary 
because on short and intermediate time scales the perturbed system is almost isolated 
from the environment. Thus, one can experimentally study if an isolated system 
equilibrates after a sharp perturbation and in this case whether it thermalizes or retains 
memory of the initial conditions. Experimental studies of the non-equilibrium dynamics 
in one-dimensional ultra-cold Bose gases have given, up to the moment, contradictory 
results [HI E] 

From a theoretical point of view, these questions have been addressed using different 
methods and points of view. The results from classical mechanics can not be translated 
directly into quantum mechanics. On one hand the concept of quantum integrability is 
not well defined, though it is generally considered as synonym for exact solvability [71 E]. 
On the other hand the unitary time evolution of quantum states leaves no room 
for a dynamical regime similar to deterministic chaos. Actually, the name quantum 
chaos stands for the different type of signatures that certain quantum systems exhibit 
depending on the large time scale behavior of their classical analogues jH |T0]. These 
signatures appear on the statistical behavior of eigenenergies and eigenfunctions and is 
not clear at all how they can influence the evolution of the system on large-time scales. 

The generally accepted assumption that integrable systems do not thermalize is 
corroborated in several models |1H [T2| [T3| [HI [15], where the non-equilibrium regime 
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extends over large time scales or where a long-time steady state with non standard 
thermal properties is reached. The approach to thermal equilibrium of generic systems 
has been studied by several authors. It has been proven that almost any system 
in interaction with a large heat bath will equilibrate and thermalize [161 (HI- Fo'^ 
isolated systems it has been shown that "typical" Hamiltonian and observables will 
be in thermal equilibrium for most times [T8l [19] . but it seems a very difficult task 
to decide whether a specific Hamiltonian or observable belongs to this class or not. 
The thermalization of specific ID and 2D fermionic and bosonic systems has also been 
studied [20| |2T| [22| |23| |2^ . Surprisingly, thermalization was not obtained in all the 
cases. Several reasons have been reported to explain the lack of thermalization in these 
systems, but it seems that the so-called eigenstate thermalization hypothesis (ETH) 
plays a fundamental role |23l |24] . The ETH states that thermalization occurs at the 
level of individual eigenstates [251 126] . whenever they satisfy Berry's conjecture on 
chaotic eigenfunctions [27j, i.e., whenever they behave as (quasi) random superpositions 
of the basis states. For this and other reasons, the role played in the viability of 
thermalization by quantum chaos in general, and by the properties of chaotic wave 
functions in particular, has began to receive some attention [281 |29] . It seems that 
the number of principal components (NPC) or inverse participation ratio (IPR), which 
keeps track of the progressive eigenstate delocalization through the integrability to chaos 
transition, might be directly related to the deviations of the steady expectation values 
from the corresponding statistical values [301 EI] 

In order to get a deeper understanding of the role played by quantum chaos in 
these processes we study the thermalization of isolated fermionic systems modeled by the 
simplest two-body random ensemble [32]. Two-body random ensembles are paradigmatic 
models to study quantum chaos and specially the dynamical transition from integrability 
to chaos. They have also been used in the past, together with some related models, 
such as nuclear shell model and interacting spin models, to perform different studies 
on thermalization. The thermalization criteria were based on the equivalence between 
different definitions of entropy [33l [Ml [35] and temperature [36]; represent ability of 
occupancies by Fermi-Dirac distribution [371 l38l [39] (Bose-Einstein distribution for 
bosons); and calculation of expectation values using the canonical distribution [38] . 
A brief review of some of these studies can be found in Refs. [lOl [H]. However, in 
the present work, as well as in most recent papers the focus is put on the Ergodicity 
principle [T8l [T9] which is the cornerstone for thermalization, and clearly more precise 
and general that the aforementioned criteria. We study the relevance of several factors 
in the thermalization process like ETH, the dimension of the Hilbert space, the structure 
and energy of the initial state and its proximity to the ground state. We also analyze 
the importance of the degree of chaos as measured by the different chaos markers. 

The rest of the paper is organized as follows. Sec. [2] briefly introduces the 
embedded Gaussian ensembles of random matrices generated by two-body interactions, 
with emphasis in the embedded Gaussian orthogonal ensemble, and some of the well 
established main features of the order to chaos transition in this ensemble. However, 
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the criteria for determining the transition points as the two-body interaction strength 
is increased, are somewhat different from those used in the past. Sec. [3] deals with the 
thermahzation of embedded Gaussian orthogonal ensembles. Here we give all the new 
results of the paper, where the transition points of Sec 12.21 play a central role. After 
introducing some basic definitions in Sec. 13.11 Sec. 13.21 reports the main numerical 
results linking thermalization with the type of spectral fluctuations, the delocalization 
of the wave functions, the structure and the proximity of the initial state to the ground 
state, or the dimension of the Hilbert space. In Sec. IHwe gather together some analytical 
results that establish a connection between thermalization, relative to a given observable 
O, and the value of the NPC for the transition strengths originated at a given eigenstate, 
or between thermalization and the dimension of the whole Hilbert space, the dimension 
of the microcanonical energy shell and the correlations generated by O. Finally, Sec. [5] 
contains the conclusions. 

2. EGOE(l+2) model: order-to-chaos transitions 
2.1. EG0E(l+2) model 

As previously mentioned, we try to analyze the relationship between the order-to- 
chaos transition and the thermalization of fermionic systems, which can be modeled 
by an appropriate two-body random matrix ensemble. Here we briefly introduce these 
ensembles and their relation with quantum chaos. Recent and comprehensive reviews 
can be found in references |40[ W\\ . 

There is a clear relationship between the energy level fluctuation properties of a 
quantum system and the large time scale behavior of its classical analogue. The spectral 
fluctuations of a quantum system whose classical analogue is fully integrable are well 
described by Poisson statistics, i.e., the spacings between successive energy levels are 
not correlated |42]. According to Bohigas, Gianoni and Schmit |43], the fluctuation 
properties of generic quantum systems, which in the classical limit are ergodic, coincide 
with those of the Gaussian ensembles (GE) of random matrices [44J. This statement, 
initially supported by many experimental data and numerical calculations, has been 
finally proven for quantum systems with few degrees of freedom, where the semiclassical 
approximation is valid [15]. The large time scale behavior of the classical analogue also 
determines the properties of the wave functions. Extensive reviews of later developments 
can be found in [T0| W\\ and references there in. 

Quantum many-body systems like complex atoms and atomic nuclei are usually 
considered to be chaotic if their spectral fluctuations are those of the Gaussian 
orthogonal ensemble (GOE), which is the appropriate GE for systems with time- reversal 
invariance and rotational symmetry. However, real quantum systems are usually well 
described by real or effective one- plus two-body interactions in the mean-field basis, 
whilst GE represent systems with multi-body interactions. The embedded Gaussian 
ensembles (EGE) of random matrices were introduced to tackle this problem, and to 



Thermalization in the Two-Body Random Ensemble 



5 



provide a more realistic picture of many-body quantum systems. Moreover, in the 
present context EGE are interesting because random interactions can illustrate the 
effects on thermalization caused by generic interactions, which lead the system from 
integrability to chaos. 

The EGE(l+2) ensembles consider m fermions or bosons distributed in n single- 
particle states \k) , k = 1,2, ■ ■ ■ ,n, interacting via the following Hamiltonian 



where the single-particle energies and the two-body matrix elements (properly 
symmetrized or antisymmetrized) {pq\ V \ kl) behave as independent Gaussian random 
variables. In this expression A gives the strength of the two-body interaction, and a|, and 
Ofc create and destroy a fermion (or a boson) in the kth single-particle state. In fermionic 
systems only the strict inequalities k < I and p < q are valid. It has been numerically 
shown that for EGE(2) ensembles (e^ = 0) the spectral fluctuations do agree with those 
of GE, provided that energy scale is redefined appropriately. Introducing the notation 
(•) = ((•)) /d, where ((•)) stands for the trace operation, and d is the dimension of 

1 /9 

the Hilbert space, the centroid and the energy span must be {E) = and {E^) ' = 1, 
respectively. [35] 

2.2. The order-to- chaos transition in EG0E(l+2) 

One of the most significant aspects of EG0E(l+2) is that as A increases, starting from 
A = 0, the system undergoes a transition from a regular to a chaotic regime that 
deeply affects the state density, level fluctuations, and wave functions. This change of 
dynamical regime is characterized by three chaos markers |40l |4T1 146| 147] . There is a first 
marker Ac that signals the transition from Poisson to GOE spectral fluctuations. This 
transition occurs when the interaction strength A is of the order of the spacing between 
the basis states that are directly coupled by the residual two-body interaction, a result 
that came out of nuclear structure calculations by Aberg [IHl HQ]. More developments 
in determining the Ac marker are given for example in [501 El]- An important outcome 
of Aberg criterion is that Ac oc l/(m?n) for EG0E(l+2), which is well verified in |51] 
and in the results presented below. 

As A increases further from Ac, the structure of the eigenstates undergoes a deep 
transformation. First the strength functions change from Breit-Wigner to Gaussian at 
a transition point denoted by Xp- Beyond \p we find a third chaos marker At which 
defines the center of a region where different definitions of thermodynamic variables, 
such as entropy, temperature, specific heat, etc., give the same results, as it occurs for 
infinite systems. As far as the statistical entropy 5'*'^^'' and the Shanon entropy 5**"-'^ are 
concerned we can understand the meaning of \t as follows. The former is proportional 
to the logarithm of the state density, which in our case has essentially Gaussian form, 
though its mean and variance depend on A. At A = 0, where the eigenstates are fully 
localized in the mean-field basis, information entropy is 5*™-^ = 0, whilst for sufficiently 
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large values of the interaction strength, the eigenstates are quite similar to those of 
GOE. In our case this occurs for A > 1 and then S**"-^ ~ log(0.48(i), except perhaps near 
the spectrum edges. Changing from A = to A = 1 the wave functions become more 
and more delocalized in the mean-field basis, being this process faster in the middle 
than in the spectrum edges, and Xt signals the precise value of the residual interaction 
strength in which the energy dependence of S'*'*^'" and 5**"-'^ is quite similar. Thus, if we 
define an appropriate "distance" between the two entropies, it should have a minimum 
at At. Another property is that At signals the duality point between h{l) and V(2) basis, 
i.e., the point where the eigenstates become equally delocalized in the two basis |17t |50]. 
Beyond At the eigenstates become quickly similar to those of GOE systems, for which 
they are essentially Gaussian random superposition of the basis states. This property 
resembles Berry's conjecture about the ergodic structure of chaotic wave functions in 
phase space [27j. Nevertheless, a word of caution should be included here. Since S*™-^ 
is basis dependent, for A S> Aj it would be more appropriate to define S**"-^ in the V{2) 
basis. 

For our purposes the markers Ac and At are the most relevant. We detect the 
change in the spectral fluctuations, and thus the position of Ac, using the nearest- 
neighbor spacing distribution, denoted P{s). The spacings for generic integrable 
systems obey the Poisson distribution, i.e, P{s) = exp (— s), while the Wigner surmise, 
P{s) = |sexp (— fs^), provides a very good approximation for GOE-like systems |44] . 
The Brody distribution [35], 

P{s,u) = A^{u + l)s'^exp{-A^s^+'), (2) 

where u is usually called the Brody parameter and A^^ is a normalization constant, it 
is used to asses how close the fluctuations are to the Poisson limit, corresponding to 
w = 0, or to the Wigner surmise with u = 1 [52]. Although Eq. ([2]) is only a heuristic 
formula, and the Brody parameter has no definite meaning for Hamiltonian systems, it 
has been employed in a variety of studies since its introduction. Very recently a physical 
foundation of u has been found by Sakhr and Nieminen in the context of self-similar 
fractals j53]. As the interaction strength in Eq. ([1]) increases from A = to a sufficiently 
large value, the Brody parameter changes from w = to u; ~ 1. In what follows, the 
position of the first chaos marker Ac is fixed by the condition w(A) = 1/2. Similar 
conditions can be found in the literature [STl 

To locate the value of the third marker At we consider the values of three different 
entropies: 

• Thermodynamic entropy, = log p{E), where p{E) is the density of states. 

• Information entropy in the mean-field basis S^^^ = — Ylk=i ^fc T where 
the coefficients cf are the eigenstates components in the mean-field basis. 

• Single-particle entropy = - X] {{ni)'^ log (^{ni)'^^ + (^l - {rii)^^ log (^1 - {rii)^^}, 
with (rii)^ the occupancy of the ith single-particle state at energy E. 
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Entropies formulas, valid for A > Xp, are given in [3U [55] for EG0E(l+2) 
ensembles. If we define the average "distance" between the three entropies as 



A,(A) 



+ {R'i-R'l^'''Y dE\ // R'l^^'^dE, 

J —oo 



(3) 



where R'^ = exp {5*^ — S!j^^^}, the value of corresponds to the minimum of A,, because 
this ensures that the values of the different entropies will be very close to each other. 




(TOOl 0.01 0.1 1 



Figure 1. Average values uj{X) (dots) and As(A) (diamonds) for A G [0, 1], calculated 
using 60 members of a EG0E(l+2) with to = 6 and n — 16. The vertical dashed lines 
indicate the position of Ac and At . 



Fig. [T] displays the ensemble averages uj{X) and As{X) along A G [0,1], for a 60 
member EG0E(l+2) with m = 6 and n = 16. In order to enlarge the region where 
the order-to-chaos transition takes place we use a logarithmic horizontal axis. The 
two vertical dashed lines indicate the respective positions of Ac and A^. It can be seen 
that Ac — 0.013 and A^ ~ 0.13, values that are consistent with theoretical estimates 
(see [10| WT\ HH]). For A ^ 0.03 we obtain u ~ 0.96, which is very close to the actual 
GOE result. However, the st ructur e of these states is still very different to tho se of a 
GOE system because 5^ ^goeI'^- O^^ly ^^en A ~ 1 one finds that ~ S'^^j^, 
signaling that the dynamics has become fully chaotic. 

3. Thermalization definitions and numerical results 
3.1. Basic definitions 

We study the thermalization properties of finite fermionic systems with time-reversal 
and rotational invariance. To be precise we consider m fermions distributed in n 
independent particle states, interacting via the Hamiltonian ([1]), where Sk and {pq\ V \ kl) 
are independent real Gaussian random variables. These systems are usually called 
EG0E(l+2) ensembles. In this work we use ^ = k, (e^ — ky = 1/2, {pq\ V \kl) = 0, 
\{pq\ V \kl)f = 1 + 5(pq),(fe«)! the energy scale is such that (E) = and {E"^)^^"^ = 1, 
regardless of the value of the interaction strength A. The results presented below have 
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m55 5 5 566 6 6 6 
n 12 13 14 15 16 12 13 14 15 16 
d 792 1287 2002 3003 5005 924 1716 3003 5005 8008 

Table 1. Matrix dimension for EGOE systems with 171 — 5,6 and n = 12 — 16 



been obtained by fully diagonalizing 60 member EG0E(l+2) systems with m = 5, 6, 
n = 12 — 16 and A G [0, 1]. The corresponding dimensions are given in Table [TJ 
Let l^'(O)) be the initial state of the system, that we decompose as 

\^m = j2c,\E,), (4) 

in the eigenstate basis {|£'^),/i= l,2,...d} of the Hamiltonian ([T]). Then, given a 
certain observable O we define: 

• The instantaneous value 0(t) = O |^(t)) = {{Op^{t))), where p^{t) = 

\m)m)\- 

• The time average (0(t))y = (2T)-i 0(r)dr. When T > 1, {0{t))^ 



{^)eq = Zl/. IC'mI Do{Ei_,), where Do{E^j) = {Ef,\ O \Ef,). The steady state average 
can also be written as (C)^^ = {{Opeq)), with peq = \Cf,f \Ef,) {E^\. 

• The statistical average {O)^^^^ = {{Opstat)), where p^tat is the density operator 
corresponding to an appropriate statistical ensemble. 

We say that the system thermalizes if for any relevant observable O and almost 
any state it is satisfied that 

m))r ^ {0).,at (5) 

In practice, to asses whether Eq. ([5]) do approximately holds for a specific 
observable, one can use the relative error 



So 









(O) • 

\^ / Stat 

which is well suited to compare the degree of thermalization of different systems, or the 
thermalization of a single system relative to distinct observables. 

Let us consider that the system is prepared in a non-equilibrium state |^1'(0)) 
for which the energy is essentially constant. Hamiltonian eigenvectors are excluded 
as they are stationary states, but we can assume that (\E'(0)| |\E'(0)) = Eq, and 
[(\I^(0)| - Eq |\Ef(0))]^^^ < AE, with AE sufficiently small compared to the energy 
spectrum span, but large enough to contain many energy eigenstates. In such a case 
the microcanonical ensemble is the preferred statistical ensemble. Denoting by W 
the corresponding energy shell, i.e., W = {\E^);E^ G [Eq — AE,Eq + AE]}, the 
microcanonical density operator is 
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Pmc = -7,^ \E^) {Ef, 



where d' is the dimension of the subspace W and the symbol ^ means that the sum 



is restricted to eigenstates belonging to W. Thus, the corresponding microcanonical 
average is given by 



To examine whether the energy plays a role in the thermalization of the system 
we shall use three different energy intervals with A^; = 0.1 and Eq = 0.0, 1.0 and 1.7. 
We denote the corresponding energy shells by Wi, W2 and W3 respectively. Moreover, 
to see how the initial conditions affect the process we let the system evolve from three 
different types of initial states, defined as: 

• |^'(^''(0)) oc Pw \ ko), where is the projector onto W, and |A;o) is the mean-field 
state with energy e = Eq. 

• |\E'(^)(0)) oc Xl^C*^ \Efj.), with the coefficients Gaussian random variables with 
mean zero and variance equal to one, i.e., = Q{0, 1). 

. |*(3)(o)) oc E^C'JE^), where the expansion coefficients are = 



The states \I'(^)(0) and \E'('^)(0) are random superpositions of the eigenstates 
belonging to W, but due to the Gaussian factor the distribution of the coefficients is 
wider for \E'(^)(0). As we shall see below the distribution of the state amplitudes inside 
W is one of the factors that may affect the thermalization process. 

Before we turn to the main results, let us introduce the four types of generic 
observables for which we study the validity of Eq. 

• diagonal one-body operators Od(l) = Ylk ^fc'^I'^fc) 

• one-body operators 0{1) = ^kidkiCi\ai, 

• two-body operators C(2) = Y.k<i,p<q^kipqa\a\aqap, 

• strength function operators Osj = 0^{1)0{1), 

where the parameters 6k, O^i and Okipq are taken as random variables. 
3.2. Numerical results 

Let us consider the ensemble averages Ao^(i), ^o{i), ^o{2), and Aq^j for a 60 member 
EGOE with {m,n) = (6,16). Fig. |2] displays their evolution with the strength A 
when the system is prepared at t = in a state ^(^)(0) with energy ii^ ~ i?o = ( 
\E'^^''(0) G Wi). Because of the large differences between the relative errors of these 
operators and to properly visualize their evolution in the short interval Ac < A < A^, 
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Figure 2. Evolution with the interaction strength A of the averages Ao^(i) (squares), 
^o{i) (diamonds), j!^o{2) (triangles), and Aq^j, (dots), given in percent, for a 60 member 
EG0E(l+2) with (m,n) ^ (6, 16) initially prepared in a state ^'(^^(O) G Wi. 



a log-log scale is used. As for Fig. [T] the two vertical lines give the positions of Ac 
and At. In all the cases Aq becomes smaller as the interaction strength increases up to 
A ~ At. It is very important to realize that the transition from Poisson to GOE spectral 
fluctuations, which is considered the most relevant signature of quantum chaos, occurs 
at A ~ Ac and does not modify this trend. On the contrary, for A > At the relative 
errors either remain essentially constant or the decreasing rate is much smaller. Recall 
that At defines a region where the three entropies S^^'^^, S^^^ and S"^ take essentially 
the same values, and signals the point at which the wave functions start to become very 
delocalized in the mean-field basis. Beyond At, Ao becomes clearly smaller than one 
percent only for two operators, namely and Ogj- Their errors are Ao^(i) ~ 0.5% 

and /S.Osf ~ 0.1%, respectively. Thus, as long as the the system is prepared in an 
initial state \E''^^)(0) G Wi and A > At, Eq. ([5]) approximately holds for the observables 
(9d(l) and Osf and then we can assert that the system thermalizes relative to these 
observables. This is not the case of the observables 0{1) and 0{2). It is worth noting 
that perhaps the main difference between 0(1), 0{2) in one hand, and 0^(1), Ogf in 
the other, is that the latter have meaningful smoothed form for large A, as given by 
spectral distribution methods [ITl [56| |57]. 

The results corresponding to other choices of ^E'(O) G Wi are shown in Fig. El For 
simplicity we have represented only the results for Osf and 0{1). Curves in black, red 
and green correspond to \I'*^^)(0), \E'*^^^(0) and \I''^^^(0), respectively. In all cases the initial 
state belongs to the energy shell Wi. We see that the choice of the initial conditions 
does not affect the main trend: Ao(i), and Aq^^ diminish progressively as the strength 
A is increased, and once A > At their values remain essentially constant. However, the 
precise values are quite different. When A > At, the initial states ^^^^^(0) and \I''^^)(0) 
give rise to very similar results while the error corresponding to \Ef*^^^(0) is clearly larger. 

To shed some light on the mechanisms leading to these results lets us note that the 
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Figure 3. Values of Ao(i) (diamonds) and Ao^j,(dots), expressed in percent, as 
function of A in a 60 member EG0E(l+2) with (m, n) — (6, 16), and initial conditions 
given by ^(^^(O) (black), vIf(2)(o) (red), and vI/(3)(0) (green). In all cases the initial 
state belongs to the energy shell Wi. 
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Figure 4. Variation with the interaction strength A of the standard deviation <^\c'^\^ 
of the initial state Vl'(O) components in the energy eigenbasis. The values of o'|p^|2 have 
been multiplied by the dimension d' of the energy shell. The results for *(i)(0) e Wi 
(black), ^'(2)(0) e Wi (red), and ¥^^0) e Wi (green) are plotted 



difference 

6o = (O),, - (O)^, = Y.' D^{E,) , (7) 

becomes very small when: the distribution of the coefficients is nearly flat, i.e, 
~ 1/'^') ^he matrix elements Do[E^) almost do not fluctuate inside the energy 
shell W , and therefore Do{E^) ^ {^)mc- This condition is known in the literature as 
the eigenstate thermalization hypothesis (ETH), which has been conjectured to hold in 
chaotic quantum systems [221 ES]- The standard deviation of the initial state components 
in the energy eigenbasis, cr|c^|2, is plotted in Fig. HI When A is very small o'^Q^^i takes 
very different values in the three cases, being quite larger for ^^^-'(0). However, in this 
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case cr|^ |2 undergoes a sharp decreasing up to A ^ Aj, where its values are very similar 
to those^of ^(2)(o) and half those of ^(3)(o). 




Figure 5. Evolution witli A of ttie matrix elements {E^\0 \Efj,) fluctuations in the 
central energy shell Wi, measured by the ratio Kq = cr™'^/ (O)^^ (see text). Results 
for 0^(1) (squares), 0{l) (diamonds), 0{2) (triangles) and Osf (circles). 



Defining the standard deviation of the expectation values Do{E^) inside an energy 
shell W as 



1 1/2 

(8) 



the ratio Hq = c""™^/ {^)mc provides us a measure of the fluctuations of these matrix 
elements. Fig. [5] plots Ko(A) for a randomly selected member of the ensemble inside 
Wi. Throughout the interval A = to A = 1 the ratio Ko is reduced sharply by a factor 
ranging between 5 and 40, but even then we get > 1 for operators like 0{1) and 
0{2). Remarkably, < 0.1 for Od{^) or Ogf, meaning that the fluctuations of their 
expectation values are small compared to (O)^^ in the energy shell Wi. 

The extraordinary resemblance of Fig. [5] with Figs. [2] and [3] shows that 
the fluctuations of the eigenstate expectation values determine for which observables 
thermalization occurs. In other words, the thermalization of the system is correlated 
with the degree of compliance of the ETH. Moreover, given a certain observable the 
exact degree of thermalization depends on the width of the initial state: the wider is 
the initial state the worse thermalizes the system, a fact that is in agreement with some 
previous claims |58) . 

Thermalization may be affected by other factors, such as the proximity of £"0 to 
the edges of the spectrum, and in particular to the ground-state energy. In order to 
explore this question we have compared the behavior of Aq in three energy shells Wi, 
W2 and W3, of width A^; = 0.1 and Eq = 0.0, 1.0 and 1.7, respectively. Although the 
energies of the largest eigenstates oscillate between 2.5 and 3 depending on the value of 
A, it is very difficult to get closer to the spectrum end. The reason is that outside the 
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Figure 6. Behavior of the ensemble averages Ao(i) (diamonds) and Ao^j (dots) as 
A varies, for three different initial state energies, Eq = (black), Eq = I (red), and 
En ^1.7 (green). 



central interval (—2,2) the state density is so scarce that it is impossible to obtain a 
narrow energy shell which contains a large number of states d', but small compared to 
the dimension d of the whole space. We plot in Fig. |6]the variation with A of the two 
averages Aq^^ and Ao(i), calculated for an EGOE system with {171,11) = (6, 16), which 
is initially prepared in a state of type ^*^^)(0). One sees immediately that Aq behaves 
very differently for these operators. Once the thermalization region (A > A^) has been 
reached, the relative error Aq^^ grows as the energy Eq approaches the end of the energy 
spectrum. By contrast, the evolution of Ao(i) is more irregular, but still this error is 
clearly larger for Eq f» 1.7 and A ~ 1. Therefore, the proximity of the initial state 
energy to the spectrum edges inhibits the thermalization of the system. 



1000 




0,001 001 0.1 



Figure 7. Variation with A of ^0(1) (diamonds), and Aq^^, (dots) in five ensembles 
with dimensions d = 924 (maroon), 1716 (blue), 3003 (green), 5005 (red), 8008 (black), 
and initial conditions given by ^'(i)(0) e Wi. 



Before closing this section we shall briefly consider whether the dimension of the 



Thermalization in the Two-Body Random Ensemble 



14 



(m, n) Hilbert space hinders or enhances thermahzation. For illustration we consider 
5 systems with (m, n) = (6, 12), (6, 13), (6, 14), (6, 15) and (6, 16), with dimensions (see 
Table ^ d = 924, 1716, 3003, 5005 and 8008, respectively. Fig. [7] shows the evolution 
of Ao for Osf and 0{1) as the dynamical regime changes from regularity at A = to 
full chaos for A ~ 1. The growth rate of d is very small, but it is enough to induce a 
gentle decrease in A^, while maintaining the same qualitative trend. Similar results are 
obtained for the other two observables Od(l) and 0{2). If we also take into account 
the analytical results of the next section we can conclude that the system will always 
thermalize in the thermodynamic limit (the difference between microcanonical and 
diagonal predictions decreases with system size), but there may be relevant differences 
for finite dimensions, which is important since many physical systems of interest are 
mesoscopic. 



4. Analytical results for thermalization 

Our goal now is to relate the relative error A^ with the fluctuation properties of the 
eigenstates and with the correlations generated by the observable O. To this end let us 
focus on the quantity 

50 = (O),, - = mpeq - Pmc))) = {{OAp)) , (9) 

with Ap = Peg — pmc- It can be considered as a random variable because we model the 
physical system by means of an appropriate ensemble and thus the external products 
\E^) {E^\ appearing in Ap change from one member of the Hamiltonian ensemble to 
another. Moreover, being interested in "typical" properties we can introduce a fictitious 
ensemble of initial states, |^E'(0)) = ^^C'/^l-E^), with similar energies. Assuming that 
they are uniformly distributed in the unit sphere in W, the fluctuation properties of 
the Cfj^ coefficients obey the Porter-Thomas (P-T) distribution [35] whenever d' is large 
enough. Thus 

where ^"w^ means averaging over the fictitious ensemble of initial states. Therefore 



Ap = implying that Sq = Sq =0. Then, Sq is the variance of ([9]) and we may 
agree on defining the "typical" value of A^ as 

1/2 



^typ 



Si 



(Of 

\ I rr, 



Introducing the shifted operator Q = (9 — a/, with a G M, it is easy to see that 
Sq = ^Q- Moreover, a straightforward calculation gives 



^ =Y^ DQ[E,)DQ[E;){\C,f - ^ (l^'^l'-^)' (^2) 
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and taking into account the fluctuation properties of tlie coefficients we obtain 
~~ 2 



(13) 



Introducing this resuh in (fTTj) . the value of A*^^ is given by the expression 

1/2 



(14) 



which can be further simplified because the averaged forms D'q^E^) are smooth functions 
of the energy that do not change substantially inside the energy window W . Thus, we 
arrive to 

n 1/2 



2 DUE,) 



d' 



DoiEo] 



(15) 



Since a is a free parameter we can choose a = Do{Eo), which gives rise to 



^typ 



2 Dl{Eo) - (^DoiEo))' 



d' 



DoiEo) 



1/2 



cr'oiEo) 



d' 



DoiEo 



1/2 



(16) 



For strength function operators, like Ogf = O^O, the value of A*^^ can be related to 
the NPC (also called IPR) in transition strengths originating from the central eigenstate 
|£^o)- Before we proceed it seems suitable to define the usual NPC for eigenstates 
and its extension for transition strengths. The former is a measure of the eigenstate 
complexity in the mean-field basis of Slater determinants \K) with K = ki, k2, ■ ■ ■ , km 
and ki = 1,2, ... ,n. Expanding the eigenstates in this basis, i.e., 



\E,) = J2 UKiE,) \K) , 



(17) 



K=l 



where the amplitudes UxiEfj) satisfy that Y1'k=i WKiE^)f' = 1, the number of principal 
components is defined as 



NPC(E) 



K 



(18) 



Note that NPC, as a function of the amplitudes U, attains its absolute minimum 
NPCmin = 1 whenever the wave function is localized in a single basis state, and its 
absolute maximum NPC,„a^x = d when the wave function is completely delocalized and 



all the amplitudes are equal to lUxiE^ 



1/d. Generally speaking NPC gives the 



effective number of basis states that build up Hamiltonian eigenstates. Thus, it is small 
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for localized states, while for chaotic states, which are quite delocalized in the mean- 
field basis, it takes values comparable but somewhat smaller than d because system 
symmetries and orthogonality prevent reaching its maximum. For instance, the average 
value for GOE is NPC''°'' = d/3. 

Since the analysis of eigenvector amplitudes is largely equivalent to dealing with 
transition strengths we can extend the previous discussion. Using the standard notation 
of spectral distribution methods, let us introduce the so called locally normahzed 
strength R, and the normalized stregth TZ generated by the action of the operator 
O on a certain eigenstate as 



R{Eu, E^) 



(EJOIE., 



\{E.\0\E^) 



(19) 

n{E,,E^)= [Mo{E,)]-'\{E,\0\E^)\\ (20) 

where Mq{E^) = (E'^l O^O \ Eij) is the total strength sum. In our notation Mo{E^) = 
Do^j:{Ee^). Note that the initial and final spaces connected by O do not need to be 
the same, but for simplicity we shall consider here a single space. If we normalize the 
state vector O and expand it in the Hamiltonian eigenbasis 

d d 

[MoiE^)]-' O \E;) = a \E.) = [^o(£^m)]"' {EA O \E;) \E,) , (21) 

v=l v=l 

then the NPC for the strength distribution generated by O on the eigenstate is 
defined as 



NPC(£;^) 



Y,(n{E^,E,))' 



(22) 



Contrary to usual NPC in wavefunctions it only depends on the operator and the initial 
state, but it is basis independent. The NPC for strength distributions gives the effective 
number of eigenstates over which the strength generated by the action of the operator O 
on a given eigenstate is spread. For operators generating collective states (then H and 
O are highly correlated) NPC should be small, while for operators generating chaotic 
states it should take large values. 

Before turning back to the calculation of A*^^ we give an expression for the 
smoothed NPC^;^ that will be used below. If we assume that the fluctuations in 
the locally renormalized strengths R{Ei,,Efj_) follow the P-T distribution then ^ = 1, 
= 3, and 



NPC, 



Y,{nE.,E,)y 



3E 



{E^OlEoY 



(Mo(£;o))" 



3j2{nE.,E,)y 



(23) 
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An EG0E(l+2) formula for NPC in transitions strengths and A not too far from 
Xt is given in Ref. |61]. The typical error A*^^ can be written as 



/vtyp 



2 DljEo) - [D^^.iEo))' 



1/2 



2M^{Eo)-(Mo{Eo) 



d' 



Mo{Eo) 



1/2 



where 



M§{Eo) = {Eo\0^0\Eor 



^ {E,\ O \Eo)' + J2 O \Eof {E,\ O lEoY 



(25) 



J2 {E,\0\Eof-(^{E,\0\Eof^ +(Mo(E, 



The step from the first to the second line in the previous expression follows from the 
independence of the strengths due to the P-T assumption, and the third line is obtained 
by adding and subtracting terms with E^ = E^. Using again that E? = ?>, this result 
can be simplified as 



M2(Eo) = 2 {{Eu\ O |i?o)')' + (Mo(i?o 



(26) 



This expression was reported for the first time in [59] (see also [IQ]), and its 
predictions have been compared with shell model results in [10| |60] for the width of 
the strength sums fluctuations. Finally, inserting (126|) into (123|) . one gets 



\typ 



4J:,{{Eu\0\Eoy 



d' 



Mo{Eo 



1/2 



n 1/2 



3d' 



(npcS™) 



-1/2 



(27) 



which establish a connection between the thermalization of the system, relative to an 
observable Osf = 0^{1)0{1), and the value of the NPC for the transition strengths 
generated by 0{1) acting on the eigenstate with energy Eq. An important outcome is 
that for chaotic systems the NPC is expected to be large and hence these system will 
thermalize, while for regular systems NPC has to be small and thus thermalization will 
be hindered. 

It is also possible to derive an expression for generic operators when X > Xf. In this 
case the eigenstates become more and more delocalized in the mean-field basis. Finally 
we reach a regime were they behave as (quasi) random combinations of the basis states, 
the squared amplitudes {Uxf abide the P-T distribution and 5**"-'^ = Sqq^ (i?*"-^ = 1). 
In this regime the expression (fT6l) . valid for a generic observable, and (127|) for Ogf, take 
on a very simple form. Using Eq. ( ITTl) to expand |£^o) in the mean-field basis, one gets 
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Do{E,) {K\ O \K) = ^ (i^l O \K) = {O) 



K=l 



K=l 



2 



{0)' + J2u'K-^U],\{K\0\K)f 



(28) 



K=l 



K,K'=1 



where the usual P-T formulas for the moments, |f/x|^ = ^/d and lUxf = 3/(i^, have 
been applied. 
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Figure 8. Ao/A^yP values for O = Od{l),0{l),0{2) and Osf in 60 member 
EG0E(l+2) systems with m — 5,6 and n = 12 — 16, prepared at t = in initial states 
*(0) = ^(^^(O) (black circles), *(2)(o) (red squares) and *(3)(o) ( 

green diamonds), 

with energy E ^ 0. 



Gathering these expressions together with f[T6|) one finally arrives to 

4 (02) T^/' 



^typ 



_dd' {oy 



(29) 



that links the value of A*^^ with the dimension d of the whole Hilbert space, the 
dimension d' of the microcanonical energy shell and the correlations generated by the 
operator (9, which are measured here by the ratio C, = [{O"^) / (O)^]^^^. Fig. |8] displays 
the values of Ao/A*^^ calculated for the four observables Cd(l), 0{1), 0{2) and Osf in 
EG0E(l+2) systems with m = 5, 6 and n = 12 — 16. These are prepared in three 
different initial states \E'(0) = \E'(^)(0), \E'(2)(o) and vI/(3)(o), with energy E ^ 0. It is 
clearly seen that < Aq/A*^^ < 3, confirming that A*^^ really behaves as a standard 
deviation. 
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Using propagation formulas for (O) and (O^) in terms of m and n one can obtain 
particular expressions of A*^^ for the operators considered in this paper. For instance, 
following the results given in Appendix E of jlQ] we have 



Od(l) 



4 /i+ 1 (^') 



1/2 

(30) 



dd' { m (e) 

for Od(l) in the so called dilute limit, i.e., m,n — > oo, m/n — > 0. Similar expressions 
can be derived for the other operators. 

An important outcome is the presence of the factor (drf')"^/^. Since the energy 
window must be sufficiently narrow so that the state density is constant inside, but 
wide enough to contain a large number of states, we can assume that d' = d/x, where x 
is essentially a fixed number. For instance, in the energy shell Wi we have x ~ 14. This 
result suggest that A*^^ is inversely proportional to the size of Hilbert space. However it 
should be noted that this trend may be modified by sharp variations of the correlations 
measured by ( (or {6"^) / (6*)^) and by the fiuctuations of the actual value Aq with regard 
to Ay. 
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Figure 9. Evolution of the relative errors A 



o<i(l), ^o(l), 



Ao(2), and 



calculated 



at A = 1, as the dimension of the Hilbert space increases. Their values have been 
rescaled so that Aq = 1 for d = 792. The dotted line stands for Aq = a/d, where the 
value a w 848 has been obtained by a least-squares fit to the data. 



Fig. [9] displays the evolution of Aq at A = 1 as the dimension of the Hilbert space 
increases. The values of A^ have been rescaled so that Aq = 1 for d = 792. In most cases 
the values of Aq fall very close to the Aq oc d~^ law, represented by the dotted line. The 
fiuctuations of the relative error Ao(2) are clearly larger due to larger oscillations of the 
corresponding ( values. This fact is largely consistent with our previous findings because 
one expects the correlations of two-body operators to be more complex than those of one- 
body operators. In order to determine more precisely if the d"^ is valid larger dimensions 
would be required, but this is a formidable computational task. Nevertheless, these 
results suggest that fully chaotic systems will thermalize relative to most observables in 
the thermodynamic limit. 
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5. Conclusions 

We have studied the thermahzation process in isolated fermionic systems, described by 
EG0E(l+2) ensembles. Our thermahzation criterion relies on the ergodicity principle 
for the expectation of observables, which is more general than other criteria previously 
used, hke the representability of occupancies by the Fermi-Dirac or Bose-Einstein 
distributions. Two-body random ensembles are paradigmatic models to study quantum 
chaos and the dynamical transition from integrability to chaos. It is well known that as 
the strength of the residual interaction is increased these systems undergo a order-to- 
chaos transition, characterized by several chaos markers like Ac, that signals the onset 
of GOE spectral fluctuations, and At which defines a region where different definitions 
of thermodynamic variables give essentially the same results, as it occurs for infinite 
systems. For A > At the Hamiltonian eigenstates become more and more delocalized 
in the mean-field basis. We have shown by means of exact diagonalizations that the 
onset of Wigner spectral fluctuations is not sufficient to guarantee thermalization in 
finite systems. Only if all chaos signatures are fulfilled, including the quasi complete 
delocalization of eigenstates, we find that thermalization occurs for certain types 
of observables, such as (linear combinations of) occupancies and strength function 
operators. 

As stated by many authors, the eigenstate thermalization hypothesis seems to be 
the mechanism responsible for thermalization, but in turn it holds only for certain 
observables provided that the eigenstates behave as quasi-random superpositions of 
basis states. We have also analyzed the influence of other factors, and found that the 
proximity of the initial state to the spectrum edges hinders thermalization. On the 
contrary large Hilbert space dimensions and initial states with a very small energy 
incertitude, i.e., with a very narrow distribution in the energy eigenbasis, enhance the 
thermalization process. 

Analytical expressions linking the degree of thermalization for a given observable O 
with different properties of the system have been deduced. For instance, we have found 
that the typical value of the relative error between the equilibrium and microcanonical 
averages, A*^^', is inversely proportional to the square root of the NPC (or IPR) for 
transition strengths generated by O acting on the middle of the microcanonical shell. 
An important outcome of this result is that for chaotic systems the NPC is expected to 
be large and hence these system will thermalize, while for regular systems NPC has to 
be small and thus thermalization will be hindered. 

Similarly, we have found that when the eigenstates become fully delocalized, Aq^'' 
is proportional to the square root of the correlations generated by the observable and 
inversely proportional to the square root of the dimension d of the whole Hilbert space 
times the dimension d' of the microcanonical energy shell . Since d' oc d, this result 
shows that Aq^^ oc d~^ and suggests that fully chaotic systems will thermalize relative 
to most observables in the thermodynamic limit. 

In conclusion we have presented the first study of thermalization in the two-body 
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random matrix ensembles using the ergodicity principle for the expectation values of 
observables. 
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